Spin echo decay at low magnetic fields in a nuclear spin bath 
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We investigate theoretically the spin echo signal of an electron localized in a quantum dot and 
interacting with a bath of nuclear spins. We consider the regime of very low magnetic fields (corre- 
sponding to fields as low as a militesla in realistic GaAs and InGaAs dots). We use both the exact 
numerical simulations and the analytical theory employing the effective pure dephasing Hamilto- 
nian. The comparison shows that the latter approach describes very well the spin echo decay at 
magnetic fields larger than the typical Overhauser field, and that the timescale at which this theory 
works is larger than previously expected. The numerical simulations are also done for very low 
values of electron spin splitting at which the effective Hamiltonian based theory fails quantitatively. 
Interestingly, the qualitative difference in the spin echo decay between the cases of a homonuclear 
and a heteronuclear bath (i.e. bath containing nuclear isotopes having different Zeeman energies), 
predicted previously using the effective Hamiltonian approach, is still visible at very low fields out- 
side the regime of applicability of the analytical theory. We have found that the spin echo signal for 
a homonuclear bath oscillates with a frequency corresponding to the Zeeman splitting of the single 
nuclear isotope present in the bath. The physics behind this feature is similar to that of the electron 
spin echo envelope modulation (ESEEM). While purely isotropic hyperfine interactions are present 
in our system, the tilting of the electron precession axis at low fields may explain this result. 
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I. INTRODUCTION 

The problem of the dynamics of an electron spin cou- 
pled by a hyperfine (hf) interaction to a bath of nuclear 
spins has been a focus of much theoretical attention, since 
the interaction with the nuclear bath is the most limiting 
decoherence mechanism in spin qubits based on quantum 
dots made of III-V materials. 1,2 While at high magnetic 
fields (B ^> 0.1 T in large GaAs dots when spin echo is 
considered) the dipolar interactions between the nuclear 
spins are the main source of dynamics leading to decoher- 
ence seen in spin echo (SE) experimenter— at lower fields 
the electron spin dephases and relaxes due to hf interac- 
tion alone. This process has been studied theoretically 
for free evolution or SE both by analytical methods 6 -^— 
and by exact numerical simulations i 10 ' 19 ~ ^£ Here we fo- 
cus on the theoretical description of SE decay due to 
hyperfine interaction alone^ i 11 ' 12 ! 16 ! 17 ! 21 since the SE ex- 
periment is currently the most developed measurement 
of coherence decay in electrically controlled gated GaAs 
quantum dotsrir— and there has been a recent progress 
in performing SE measurements on optically controlled 
electron spins in quantum dot o 27 ' 28 and spins of electrons 
bound to donors^ 

Most of the analytical approaches were concentrated 
on the "perturbative" regime^^ 3 - -15 ' 18 of magnetic fields 
in which the electron Zeeman energy and the total hf 
energy A = J2i Ai (where the sum is over all the nuclei 
and Ai are the individual hf couplings) fulfill _4/il<C 1. 
Only recently it has been propose d 16 ' 17 that an analytical 
calculation can be well-controlled under a much weaker 
condition of 5 = A/QVN <C 1 (N being the number 



of nuclei interacting with the central spin) , allowing the 
calculation of decoherence at much lower B fields (which 
only have to be larger than the typical Overhauser field 
due to the nuclei, which is on the order of a few mT in 
large GaAs dots). In this theory, as in earlier closely re- 
lated studies^ ' 11 ! 12 an effective pure dephasing Hamilto- 
nian describing hyperfine-mediated interactions between 
the nuclei is used. For a pure dephasing problem one 
can formulate a diagrammatic expansion technique for 
spin decoherence time evolution, and a class of diagrams 
of leading order in 1/N expansion can be resummedjii 
leading to predictions for narrowed stat o 30 ' 31 Free Induc- 
tion Decay (FID) , SE decay, and also decoherence under 
any other dynamical decouplin g 23 ' 32 " — sequence of ideal 
7r pulses driving the qubit. 

In this article we present a comparison between the SE 
decay calculated using the ring diagram theory (RDT) 
of Refs. and exact simulations of a system with 

A" = 20 nuclei. In this way we set out to clarify the limits 
of quantitative applicability of the RDT, i.e. the ranges 
of electron spin splitting and time-scales on which the 
analytical theory based on effective Hamiltonian accu- 
rately describes the SE decay. We also investigate nu- 
merically the regime of very low spin splittings, at which 
the RDT is bound to fail, and we analyze simplified ap- 
proaches which can be used to model (at least on a certain 
timescale) the SE signal in this regime. 

The exact results confirm that the RDT is quantita- 
tively accurate as long as 5 -C 1 , and that the presence of 
multiple nuclear isotopes (having different Zeeman split- 
tings) is crucial for qualitatively correct description of 
the SE decay at these low fields. While with N = 20 
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spins used in the exact calculation it is impossible to 
unequivocally discern whether it is 6 <C 1 or A/£l <C 1 
which is controlling the quantitative agreement, our sim- 
ulations strongly suggest that the qualitative (or even 
semi-quantitative) predictions of the RDT still hold even 
at lower B fields. We observe that RDT correctly pre- 
dicts many qualitative features of the dephasing process 
even for <5 ~ 1, i.e. the RDT results are valid (at least 
qualitatively, and at least at certain timescale tr) beyond 
the regime A/Vl <C 1. The striking qualitative difference 
in SE decay between a heteronuclear and a homonuclear 
system is clearly visible for <$ ~ 1 (when A/ CI > 1), and 
the timescale at which the majority of the decay (say, by 
half of the initial amplitude) occurs is reproduced by the 
RDT. This leads us to the conclusion that the smallness 
of A/O. is not necessary for correct analytical description 
of the spin echo signal. 

Another interesting feature of the low--S SE signal in a 
homonuclear system is the oscillation of this signal with 
time at a frequency given by nuclear Zeeman splitting. 
This feature is clearly visible in our results, and it should 
not be confused with the oscillations appearing in a het- 
eronuclear situation, when the frequencies are given by 
the differences of nuclear Zeeman energies of different 
isotopes, 1 ^ 7 .^ 

The paper is organized in the following way. In Sec. [TT] 
we give a brief outline of the analytical RDT theory and 
its limitations, and we describe the system of 20 nuclei to 
which an exact numerical method of evaluation of spin 
dynamics is applied. In Sec. IIIII we present the results 
of the calculations, and in Sec. |IV]we discuss them and 
compare them with simplified "box wavefunction" model 
(in which all the hf couplings are the same). In Appendix 
|A1 we provide an extended discussion of the influence of 
choice of the hf couplings on the results of the calculations 
in a system with N = 20 nuclei. 



II. THEORETICAL APPROACH 

The Hamiltonian is given by 

i i 

where f2 is the electron spin Zeeman splitting, uj a u] is the 
Zeeman energy of the i-th nuclear spin which belongs to 
nuclear species a, and the last term is the hf interaction. 
We employ J— 1/2 nuclear spins in the paper. 



A. Overview of the effective-Hamiltonian-based 
theory of ring diagram resummation 

When 17 is large enough one can perform a canoni- 
cal transformation on the full hf Hamiltonian given in 
Eq. (TT]), which removes the electron- nuclear spin flip 
terms (5 J?) in favour of a hierarchy of intra-bath 



hf-mediated interactions involving two or more nuclear 
spins &iL£ r 17 i 21 It has been argue d 16 ' 17 that as long as 
5 -C 1 and when considering the evolution up to a certain 
time-scale tr, one needs to take into account only the 
lowest-order term in the the expansion of the effective 
hamiltonian H: 

where the first two terms are renormalizations of the nu- 
clear and electron Zeeman energies, respectively, and the 
third term is the hf-mediated interaction. 

With such an effective pure dephasing Hamiltonian 
(where the only electron spin operator present is S z ) 
the diagonal elements of the reduced density matrix 
of the central spin (giving the average spin in the 
z direction) remain constant in time, and the evolu- 
tion of the off-diagonal element p^ (t) can be mapped 

onto evaluation of the bath-averaged contour-ordered 
exponent 1 12 ' 16 ' 17 i 36 ~ — This formulation of the problem 
allows one to employ some of the tools of the diagram- 
matic perturbation theory, most importantly the linked 

cluster theorem, using which one can write p^ as an 

exponent of the sum of all the connected diagrams con- 
tributing to the original expansion. 

The third term in Eq. ([2]) is a long-range interaction 
coupling all the N nuclei within the bulk of the electron's 
wavefunction with comparable strength. For such an in- 
teraction all the terms in the perturbative expansion of 

p^ (f) can be classified by their 1/N dependence. The 

leading-order terms can then be resummed and the time- 
dependence of the SE signal can be easily calculatedj 16 ' 17 
Let us stress that both the exponential resummation of 
the original expression for p^ (t) and further resumma- 
tion of a class of linked diagrams are technically feasible 
because we use an approximate effective pure dephasing 
Hamiltonian. 

The decoherence due to the bath dynamics caused by 
the dipolar interactions between the nuclear spins (spec- 
tral diffusion) can can be calculated in a similar way, 
but only the second order (in dipolar coupling) diagram 
needs to be retained in the cluster expansion in order to 
get a controlled description of SE decay ; 5 ' 6 ' 12 and a very 
good agreement with spin echo measurements of spins 
of phosphorus donors in Si was obtained^ 5 - ' 39 ' 40 Let us 
also mention that the theoretical prediction 4 - - — of the co- 
herence decay timescale of > 10 fis due to the spectral 
diffusion at high B in GaAs quantum dots was recently 
confirmed^ 6 - According to the RDT, in GaAs dots at sub- 
Tesla magnetic field the hf- mediated interactions give the 
SE decay which is an order of magnitude faster than the 
decay due to dipolar interactions) 16 ' 17 and since we are 
interested in correct description of coherence dynamics 
on timescale comparable to its characteristic decay time, 
we neglect the dipolar dynamics altogether in the low-B 
regime in which our focus is on relatively short times. 

For the relevant here case of spin echo we define the 
decoherence function W(t) = p^ (t) (assuming p^ (0) = 
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1), with which the expectation values of the transverse 
components of the central spin are given by (S x (t)) — 
\ReW{t) and (S v (t)) = -§ImW(*). 

The RDT calculation of decoherence due to two-spin 
hf-mediated interactions simply involves diagonalization 
of NxN matrix at each time step. Such a T-matrix can 
be easily calculated, as described in details in Ref. |l7|, and 
the decoherence function is given then by the following 
formula involving the eigenvalues A; (t) of the T-matrix: 



N 



w(t)=n 



i 



W) 



(3) 



Let us note that the pair correlation approximation 
(PCA) from Refs. iEM! can be extended to the hetero- 
nuclear case considered here, and it leads to approximat- 
ing the decoherence function by W / pcA = exp(— | J2itf )- 
This corresponds to taking only the lowest order term 
in the linked cluster expansio n 12 ' 17 of W(t), while in the 
RDT we resum all the terms in this expansion which are 
of the leading order in 1 /N at every order in the spin-spin 
coupling. 

The key feature of the spin echo sequence in the high- 
field < 1) regime is that the dephasing of the cen- 
tral spin due to .^-conditioned interaction from Eq. ([2]) 
with a single nuclear species (i.e. all the nuclei having 
the same Zeeman splitting) is practically completely un- 
done by the pulse sequence.- 6 -"^ However, at magnetic 
field lower that B c ps ^JQ/ AujA/ VN\g c s |/Xb (where Aw 
is the typical difference of Zeeman splittings between dis- 
tinct nuclear species and g c g is the effective g- factor of 
the electron) the hf-mediated processes between nuclei of 
different species were predicted to completely dominate 
the SE decay . 16 ' 17 This prediction has recently been con- 
firmed experimentally in a spin echo experiment in a dou- 
ble dot singlet-triplet qubit^ 6 - in which the B-dependcncc 
of the SE decay and the characteristic oscillations of the 
SE signal with frequencies oc uj a p = uj a — ujp have been 
seen. Let us note that with the parameters used in the 
calculations below the B field regime in which these het- 
cronuclear process dominate the decay corresponds to 
£l<£! r «100 in the units defined below. 



B. Approximations inherent in the ring diagram 
calculation 



There are higher-order terms in the expansion of the ef- 
fective Hamiltonian H, which correspond to virtual pro- 
cesses involving more than two S-J spin flips These 
terms also involve increasing numbers of nuclear spins, 
i.e. there are n-spin interactions in H^ n \ while the RDT 
can only deal with two-spin interactions. The exact 
derivation of these terms quickly becomes very compli- 
cated, but we do not expect such an expansion to be use- 
ful: when the first higher-order interaction in H becomes 
important, all the higher order terms become equally rel- 
evant. However, in order to illustrate the breakdown of 



the RDT we will use in the calculations the two-spin S z - 
independent interaction which appears in H^ 3 \ in which 
the coupling constants are proportional to A 3 /N 2 £l 2 . 

A more important approximation underlying RDT is 
the fact that the effective-Hamiltonian approach neglects 
the transformation of states associated with the canon- 
ical transformation (i.e. only the Hamiltonian is trans- 
formed). As a result, no decay of S z component of the 
central spin can be obtained. This also has consequences 
for the decay of the transverse spin components at low B 
fields, where the effective-Hamiltonian based theory does 
not reproduce the fast "visibility loss" process* 6 -!*^ The 
latter process, in which the transverse spin component 
decays by ~ 5 2 on a time-scale of y/~N/A can be cap- 
tured by an exact simulation or by a theory in which at 
least some effects of state transformation are retained, see 
e.g. Ref. H. A more detailed calculation of this process, 
which is associated with decay of the S z component of 
the central spin, is given in Ref. [9|, where the Generalized 
Master Equation method with the full hf Hamiltonian 
was used. 

The transformation of states corresponds to an entan- 
glement of electron spin with the nuclei (each possible 
initial bath state gets entangled on a short time-scale 
with the central spin state). Semiclassically it can be 
envisioned as a slight tilting of the electron quantization 
axis from the z direction, which happens due to hf inter- 
action with the collective spin of the nuclei. The latter 
is semiclassicaly reproduced as a random vector, and the 
tilting of the central spin's z-axis is conditioned upon the 
size and the direction of this vector. 

The visibility loss is only one example of effects which 
are absent in an effective-Hamiltonian approach. There 
are also other features, which require the treatment based 
on the full hf Hamiltonian. The analytical approaches 
of this kind have only been applied to the case of FID 
decay in a narrowed nuclear state.— ' 13 ' 14 ' 18 At times much 
longer than N/A an asymptotic 1 jt 2 coherence decay was 
obtained ) 13 ' 14 ' 18 which is not reproduced by the RDT 
calculation! 16 ' 17 However at high fields (A/Q < 1) most 
of the narrowed FID decay is described by an exponential 
e -t/T 2 obtained both by the RDT (and PCA at very high 
fields 1 - 1 -) and by the Generalized Master Equation (GME) 
approach using the full Hamiltonian in Ref. 18,, and using 
the effective Hamiltonian in Ref. 

It should be stressed that in the full-Hamiltonian the- 
ories it was argued that the convergence of the calcula- 
tion is guaranteed only when A/ SI is small, which is a 
much more restrictive requirement than i < 1 required 
in RDT theory. However, neglecting the higher-order 
multi-spin interactions in the effective Hamiltonian was 
arguedi 7 - to be a good approximation for S <C 1 only up 
to a certain time-scale r^j, so that even if the effects of 
the state transformation were unimportant, the temporal 
regime of applicability of RDT is certainly limited. On 
the other hand, the GME theories seem to give a well- 
controlled result for all times, but in a more restricted 
regime of magnetic fields. Let us also mention that tr 
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might be different for narrowed FID and SE, so that the 
conclusions on the relation of our results on SE presented 
here to the results obtained for narrowed FID with other 
theories should be drawn carefully. 

An experimentally relevant question is whether tr is 
larger than the timescale of significant coherence de- 
cay. While the estimate of tr ~ N/A given previ- 
ously is currently enough for making contact with recent 
experiments 2 ^— on spin echo in GaAs where it corre- 
sponds to about 10 /is, establishing a more precise bound 
is of large current theoretical and possibly future exper- 
imental interest. 



C. Numerical treatment of the full hyperfine 
Hamiltonian 

In order to establish more firmly the regime of applica- 
bility of the RDT one should compare its predictions with 
exact calculations starting from the full hf Hamiltonian 
from Eq. (p} . In this way it is also possible to investigate 
the regime of very low B fields (corresponding to 8> 1), 
in which no comprehensive analytical approach has been 
successful. In the exact numerical simulation the time- 
dependent Schrodringcr equation for the central spin and 
the bath is solved using the Chebyshev polynomial ex- 
pansion of the evolution operatori 20 i 22 

Although the nuclear spins in GaAs and InAs have 
J = 3/2 or 9/2, we use J = 1/2 in the paper. The main 
expected effect of such a simplification is slowing down of 
the decay (larger nuclear spins are more efficient at deco- 
hering the central spin). It allows us however to use the 
maximal number of nuclei N = 20 in an exact simulation, 
diminishing thus the possibility that differences between 
the numerics and RDT simply come from the failure of 
1 /N expansion of diagrams describing various processes 
in the nuclear bath. 

In a real quantum dot the number of nuclei interact- 
ing appreciably with the electron is commonly defined as 
N = A 2 / J2i A 2 - We identify this quantity with the finite 
N used in the simulations. This leads to 6= A? /ft. 
Furthermore, we use energy units in which J^A? = 1, 
so that S is equal simply to I/O and A = in these 
dimensionless units. The time is then measured in units 
of VN /A, which in real GaAs dot with N rj 10 6 corre- 
sponds to about 10 ns (with J =1/2 the T 2 * decay time 
of inhomogeneously broadened FI D 10 ' 41 is equal to \/8 in 
the dimensionless units used in the figures). 

Our calculations correspond to the following experi- 
mental procedure. Initially the electron spin is assumed 
to be directed along the x axis in the |+a;) state. Af- 
ter evolution for time t/2 an instantaneous rotation by 
7r about the x axis is applied to the central spin. The 
expectation value of the S x operator at the final time t 
is then the SE amplitude, plotted in all the figures as a 
function of the total pulse sequence time t. The general 



formula is given by 

(S x (t)) = ( (+x\ e ikr a x e iilr ^-e- i6r a x e- i6T \+x) ) 

(4) 

where (...) = iTrj[...] is the average over an unpolarized 
nuclear bath (with Z = 2 ), and r = t/2. 

The calculations have been performed on a model sys- 
tem of N = 20 nuclear spins coupled to the central spin 
by hf Hamiltonian from Eq. (JXJ) . The hf couplings Aj 
were drawn from a random distribution, with constraint 
Eili^?-1- We discuss the choice of such a distribu- 
tion, and also the sensitivity of the results to the specific 
choice of Ai, in Appendix [A] Here we remark that with 
N = 20 the SE signals are very similar one to another for 
most the choices of the set of Ai, and below we employ 
the set of Ai corresponding to a rather typical SE signal. 
In the hetero-nuclear case we have divided the nuclei into 
3 species with 10, 6, and 4 members in order to mimic 
the ratios of concentrations of various isotopes in GaAs 
( 75 As, 69 Ga, and 71 Ga, respectively). The Zeeman split- 
tings uj a of these 3 species were fixed at 0.02526, 0.0354, 
and 0.045, again mimicking the ratios of nuclear Zeeman 
energies in GaAs. Note that the uj a are kept fixed while 
the electron Zeeman splitting Vt is varied. Although this 
does not correspond to realistic situation when all the 
Zeeman energies are proportional to the B field, it allows 
us to more clearly separate the effects that the electron 
and nuclear Zeeman energies have on the time depen- 
dence of the SE signal. Let us note that while the ratio 
f2/cj~ 10 3 in GaAs, here we are using ft/cu^ 10 — 10 2 for 
the range of considered below. Another difference in 
comparison with the realistic GaAs dot is that the differ- 
ences of hf couplings Ay = A, — Aj (i.e. the differences 
of Knight shifts of different nuclei) are of the same order 
of magnitude or larger than the nuclear Zeeman energy 
differences wy, while in GaAs dot with N> 10 5 we have 
Aij <C LJij in the whole range of B fields for which d<g. 1. 
Thus, unless we put all the Aj equal one to another, we 
do not expect to see the w^-related oscillation in the SE 
signal in the heteronuclear system (see Ref. |l7| for details 
on why the Ay -Cwy condition is needed to obtain these 
oscillations). The existence of this oscillation has been 
confirmed by recent experiments^ and here we focus on 
other aspects of SE decay dynamics. 



III. RESULTS FOR SPIN ECHO 

First, let us briefly recount what we expect to see in 
exact numerics at not-too-low magnetic fields (5 <C 1) 
based on previous analytical work on spin echo. There 
should be a "visibility loss" initial decay of SE signal of 
amplitude ~ S 2 occurring at timescale of y/N /A (i.e. of 
order of 0(N ) in units employed here), and we expect 
a qualitative difference in the magnitude of SE decay 
between the homonuclear and heteronuclear bath. The 
open questions are: what is the timescale tr on which 



5 




100 200 300 400 500 600 

t 



FIG. 1: (Color online) SE decay for N = 20 nuclei calculated 
exactly for fi = 0.1, 1, 2.5, 5.5. The hf couplings Ai are given 
as set 3 in Table [I] The energy units are such that £\ A\ = \. 
The corresponding time units are such that T2 = y/B (for 
comparison, T 2 * w 10 ns in GaAs dots) . 



the RDT remain quantitatively accurate for 6 -C 1, and 
what happens at very low B fields at which 8 ~ 1 or even 
r5 > 1. 

In Figure Q] we are presenting the results of exact cal- 
culation for a hetero- nuclear bath with electron Zeeman 
energies fie [0.1, 5.5] (corresponding to S& [0.18, 10]). As 
expected, the exact calculation shows a very fast decay 
having S 2 magnitude for small S. In Figure [2] we com- 
pare the exact results for fi = 2.5 and 5.5 with the RDT 
calculations. These calculations were done using the S z - 
conditioned interaction from Eq. @, and also with the 
S^-independent two-spin interaction appearing in the 3rd 
order of the expansion of the effective Hamiltonian,— in 
which the coupling constants are smaller by A/ Nil fac- 
tor. For fi = 2.5 one can see that this 3rd order interaction 
is not completely irrelevant, signifying the importance of 
higher-order corrections to the effective Hamiltonian ap- 
proach. At even lower values of fi the RDT calculation 
fails to quantitatively describe the decay : not only it 
does not capture the very fast initial drop, but also at 
longer times it predicts a decay much faster than the 
one given by the exact calculation (these results are not 
shown, but the beginning of such disagreement between 
RDT and the exact signal can be seen at fi = 2.5). The os- 
cillatory character of the fi = 0.1 signal will be discussed 
later in the paper. 

On the other hand, at a slightly larger field fi = 5.5 
(when S = 0.18), the RDT calculation using only the 
lowest-order hf-mcdiated interaction is approximating 
very closely the exact result, and the higher-order cor- 
rections are irrelevant. The comparison between RDT 
and PCA in this case shows how the resummation of all 
the ring diagrams extends the timescale on which the 
analytical theory closely matches the exact calculation. 
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FIG. 2: (Color online) (a) Comparison between the exact 
(solid line) result for Q — 2.5, the result obtained using the 
RDT (red dashed line for 2nd order 5' z -conditioned interac- 
tion from Eq. J2]), blue dashed line for ^-independent two- 
spin interaction for the third order effective Hamiltonian), 
and the result obtained using the Pair Correlation Approxi- 
mation from Ref. (dot-dashed line), (b) same as (a) only 
with fi = 5.5. 



With N — 20 it is hard to say with full confidence whether 
it is the smallness of A/£l or 6=A/v~NQ which controls 
the agreement between the RDT calculation and the ex- 
act result. However, it is clearly visible that the onset 
of long-time agreement between the two calculations cor- 
relates with the suppression of the short-time visibility 
loss, which is known to be controlled by 5. This strongly 
suggests that it is indeed the smallness of S that makes 
the RDT work. Finally, these results are showing that 
the timescale on which the RDT calculation of SE sig- 
nal is quantitatively accurate visibly exceeds a value of 
\/20. which was conservatively estimated in 

Ref. GJ. 

In order to illustrate the degree to which the low-field 
SE decay is dominated by processes involving nuclei of 
different species, we have performed the calculations for 
the homonuclear bath, in which we set all the nuclear 
Zeeman energies equal to a common value u. In Fig. [3] 
we show the exact calculation for fi = 1, 2.5, and 5.5. 
Now instead of the decay reaching \ W\ <0.1 in Fig. [T]we 
see only the initial visibility loss followed by an oscillation 
with frequency given by the nuclear Zeeman splitting lu. 

Within the RDT framework, the lowest-order term in 
H which contributes to the SE decay in the homonu- 
clear bath is the S z -independent two-spin interaction 
from H( 3 \ The decay calculated by RDT with this in- 
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FIG. 3: (Color online) SE decay calculated exactly for iV = 20 
with parameters as in Fig. [T] only with all nuclear Zeeman 
energies set to a; = 0.0354. 
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FIG. 4: (Color online) Comparison of the spin echo decay in 
a heteronuclear bath (dashed lines, parameters the same as in 
Fig. [1} and a homonuclear bath (solid lines, parameters the 
same as in Fig. [2]) for SI = 2.5, 1, 0.1 (with the initial decay 
smallest for the largest SI). 



between this oscillation and the well-known electron spin 
echo envelope modulation^ (ESEEM). 

In Fig. U we show the comparison between the exact 
SE decay in a hetero- and homonuclear bath at very low 
magnetic fields (SI < 2.5), at which the RDT fails quan- 
titatively. One can see that the qualitative difference 
between the SE decay in the two cases (homo- vs het- 
eronuclear bath) is also clearly visible for Q = 1, which 
is already completely outside the domain of applicabil- 
ity of RDT. In this strongly non-perturbative regime the 
magnitude of the initial decay ceases to be equal to 5 2 , 
and instead it is smaller (with |VF|^0.7 after the initial 
decay). For the homonuclear bath we then have a small- 
amplitude oscillation about this value, while the decay 
in the heteronuclear case is practically complete on the 
considered timescale. At longer times this decay is very 
similar to the signal for f2 = 2.5, which also shows the 
robustness of qualitative predictions of RDT outside the 
regime of its quantitative applicability: within RDT one 
obtains _B-independent SE decay below a certain value of 
magnetic fields At even lower values of SI the situation 
changes: for SI = 0.1 the homo- and heteronuclear case 
signals start to look qualitatively similar. Although the 
physical picture of hf- mediated interactions is not strictly 
applicable in this regime, one could qualitatively describe 
this behavior by saying that the higher order hf-mediated 
interactions (beyond the second order one, for which the 
homo- and heteronuclear baths give qualitatively differ- 
ent SE decay) become very strong, and the difference 
between the two cases disappears. Thus, the oscillatory 
signal at SI = 0.1 has a common origin in homo- and 
hetero- nuclear cases, and the discussion of it is given in 
the next Section. Let us also note that the SI = 0.1 SE 
signal for the heteronuclear case does become negative at 
some times. 



IV. DISCUSSION AND COMPARISON WITH A 
SIMPLE MODEL 



teraction for fl = 2.5 in the homonuclear case is very 
similar to the result shown in Fig. [2^, as blue dashed line, 
and it does not show any oscillation. The same holds 
for fl = 5.5 case, where RDT predicts practically no de- 
cay in the homonuclear case, while the exact calculation 
gives the signal oscillating with the nuclear ui frequency 
around the value of 1 — 8 2 . Similar oscillations have been 
found before in numerical simulations of the multi-pulse 
dynamical decoupling protocols for an electron spin de- 
cohered by the nuclear spin bath, 23 where full hyperfine 
interaction was considered. This shows that this oscilla- 
tion is a feature following from the full S • J, interaction 
between the electron and the nuclear spins, and as such 
cannot be captured by the effective Hamiltonian theory. 
In this it is similar to the fast initial oscillation with fre- 
quency ~ SI accompanying the visibility loss. In the next 
section we discuss a similarity (also noticed in Ref. |23| ) 



The presence of the oscillation with ui frequency in a 
homonuclear system can be derived using a simplified ex- 
actly solvable model ffi 10 ' 43 in which all the hf couplings 
Ai are put equal to the same value A — 1/y/N. The 
hf Hamiltonian is then given by AS -J with J = Y). Jj 
being the operator of the total spin of the N nuclei. 
The calculation of the evolution of the system is most 
straightforward using the basis of eigenstates of J 2 and 
J z . The only pairs of states coupled by the Hamiltonian 
are then |±,j, m) and|=p, j,m± 1), with the first quan- 
tum number specifying the central spin state (± for elec- 
tron spin up/down), and we have J 2 \j, m) =j(j+l) \j, m) 
and J z \j, m) = m\j, m). Such a simplification of the 
dynamics of the system does not occur in a heteronu- 
clear case. There we have the full Hamiltonian given 
by H = VlS z + E a w a j* + £ a A Q S • J Q , and even if 
we put A a — A we must retain distinct oj a . If we then 
look at the dynamics starting with the initial state, say, 
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FIG. 5: (Color online) (a) The comparison between the ex- 
act results (solid lines) for spin echo in a homonuclear system 
(with the nuclear Zeeman energy a; = 0.0354) for electron spin 
splitting Q. — 2.5, 0.1 and the results obtained within the 
model with all the hf couplings Ai being the same (dashed 
lines), (b) The comparison of such box-wavefunction calcula- 
tions with different numbers of nuclei N for £1 = 2.5. The fast 
oscillation is an artifact of small N and the box model, and 
it disappears with increasing N. 



|+, {j a , Tn a }) (described by a set of nj pairs of quantum 
numbers j a ,77i Q for a=l. ..nj), we see that the Hamilto- 
nian couples this state to a family of states with central 
spin down and one m a increased by 1. Now, unlike in the 
homonuclear case, these states can also couple to other 
states which have the electron spin up and mp (with 
(3 ^ a) decreased by 1. Thus, while in the homonuclear 
case we only have to solve multiple two-state problems 
to obtain the system dynamics, in the heteronuclear case 
we still have to consider Hilbert spaces of dimensions as 
large as 2 Ila^a + 1) where N a is the number of nuclei 
of a species. 

The details of this calculation are given in Appendix IB1 
here we focus on the result shown in Fig. [5] In Fig. [SJi we 
present a comparison between the exact calculation with 
nonuniform Ai (solid lines) and the uniform Ai calcula- 
tion (dashed lines). The slow oscillation with frequency 
~ uj and the amplitude of the signal is reproduced by 
the "box wavefunction" calculation. The fast oscillation 
(with frequency proportional to f2) visible in the constant 
Ai calculation is the small N artifact related to artificially 
regular structure of the Hamiltonian spectrum when all 
Ai are the same. In Fig. [SJd we show the uniform Ai 
results for N = 20, 10 3 , and 10 4 . At larger N the fast 
oscillations disappear, and the shift of the oscillation fre- 
quency, while visible, saturates quickly (the results for 
N = 10 3 and 10 4 are practically the same). This allows 
us to conclude that the presence of the w-oscillation is 



not a small N effect, and that such an oscillation should 
be present in the SE signal of a central spin interacting 
with a large homonuclear bath. 

For large N the uniform Ai model corresponds to the 
situation in which a central spin S interacts with a large 
spin J. This leads us to a natural conjecture that the w- 
oscillation follows from classical dynamics of the total nu- 
clear spin, i.e. it can be recovered from classical equations 
of motion for two spins coupled by isotropic Heisenberg 
coupling and each experiencing a different Zeeman split- 
ting. Since the typical magnitude of the classical nuclear 
spin J is proportional to yN, the electron spin precession 
due to hf coupling is much faster than the nuclear spin 
precession due to interaction with the electron spin. 41 
This leads to averaging out of the electron-induced nu- 
clear precession, which leaves only the nuclear spin pre- 
cession due to the Zeeman term. One is then looking at 
a problem of electron spin dynamics due to the magnetic 
field and also due to the hf coupling with the classical 
J vector precessing with frequency to about the z axis. 
The appearance of u frequency in SE signal is then nat- 
ural. Our preliminary calculations involving averaging 
of the classical dynamics of coupled S and J spins over 
the ensemble of initial values of Jo (without making the 
adiabatic approximation used above) support this con- 
jecture. This conjecture is also in agreement with the 
physical picture suggested to explain similar oscillations 
for multi-pulse dynamical decoupling protocols^ 

Here we can comment on the f2 = 0. 1 result from Figs. [1] 
and|4j At such a low electron spin splitting the SE signal 
oscillations are present also for the heteronuclear bath. 
These oscillations are controlled mostly by the 3 nuclear 
Zeeman frequencies (which we checked by varying ui a , 
not shown in the Figures), suggesting that the physical 
picture described above, albeit with 3 classical nuclear 
spins, might be a starting point for the description of 
very low B behavior of the spin echo. 

For Q 1 one can also see a close connection between 
w-oscillation of the SE signal and the ESEEM, discussed 
previously^ in the case of a central spin interacting via 
anisotropic hf interaction with the nuclei^ in which case 
the J X S Z coupling was a source of an oscillation of the 
SE signal. In our case, if we neglect the "visibility loss" 
effect, the effective Hamiltonian for tt^> 1 is given by the 
sum of H = ttS z + uj z + AS Z J Z and the hf- mediated 
interaction from Eq. ([2"|l. for which we obtain a perfect 
recovery of coherence in SE experiment. However, the 
initial "visibility loss" corresponds to a tilt of the electron 
spin quantization axis from the original z direction by an 
angle proportional to 5, and thus the decrease of (S z ) by 
a factor ~ S 2 . If we then rotate the coordinate system 
so that the new z' direction is along the new electron 
spin quantization axis, we will obtain the terms of S zl J xl 
type in Hq defined above, thus arriving at the problem 
analogous to ESEEM due to anisotropic hf interactions. 

The "box wavefunction" approach was used previously 
to calculate \ow-B free evolution of S z component of the 
central spin^ and it was shown to agree with the exact 
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FIG. 6: (Color online) (a) Comparison of the SE signal cal- 
culated with non-uniform hf couplings (red and blue solid 
lines) and with all Ai — A (dashed lines for N = 20, circles 
for N — 40) for homonuclear and heteronuclear baths with 
N = 20 and at fi = 5.5. The thin black solid line is the RDT 
result for uniform couplings (which is N independent), (b) 
The same for Q, = 2.5. 



solution (for a system with inhomogeneous Ai) on the 



timescale of T 2 * ~ y/N / A. Above we have seen that the 
box model captures the main features of the SE signal in 
a homonuclear bath on a much longer timescale, which 
exceed not only T| but also N / A. 

Interestingly, this agreement between the "box" model 
and the exact calculation breaks down at shorter 
timescale for the heteronuclear bath. In Fig. [5] we show 
the comparison between the nonuniform Ai calculations 
and the "box" results for homo- and heteronuclear baths 
(at Q = 2.5, 5.5 and with N — 20, 40 nuclei). The RDT 
results obtained with the uniform hf couplings are also 
shown (this calculation gives the SE signal independent of 
N) . For fi = 2.5 there is an obvious difference in long-time 
behavior between the "box" calculation and the calcula- 
tion with nonuniform Ai . The "box" calculation exhibits 
a saturation of the SE signal at t ~ 400 (with the sig- 
nal staying close to this value for much longer times, not 
shown in the Figure). However, the amplitude of this 
saturation decreases with increasing N. While we have 
2(S' a; ) sat «0.46 for 7V = 20, for iV = 40 we obtain a value 
of about 0.22, suggesting that these results are affected 
by finite- N effects. The reason for which the these effects 
are spoiling the agreement between the exact and "box" 



calculation only in the heteronuclear case remains to be 
further elucidated. 

At higher f2, as in Fig.|BJi, the finite- N effect is weaker, 
and with N — 40 the RDT and the full- Hamiltonian cal- 
culation agree very well for uniform Ai. Both of these 
signals exhibit a somewhat faster decay than the exact 
result with nonuniform couplings (which is well repro- 
duced by the RDT calculation with nonuniform Ai, as 
shown in Fig. [2b). This shows that at iV = 20 the choice 
of the specific set of Ai couplings can have a visible im- 
pact on the SE signal. This is discussed in more details 
in Appendix [A"l 



V. CONCLUSIONS 

It this paper we have investigated the low magnetic 
field spin echo (SE) signal of an electron spin inter- 
acting via the hyperfine (hf) coupling with the nuclear 
bath. Intrabath dipolar interactions have been neglected. 
We have used three theoretical approaches: the exact 
numerical solution for a relatively small nuclear bath, 
the analytical theory based on resummation of ring 
diagrams ; 16 i 17 and calculations in which uniform hf cou- 
pling with all the nuclei was assumed. 

The exact numerical calculation strongly suggests that 
the ring diagram theory (RDT) describes quantitatively 
the SE decay when the electron Zeeman splitting fl is 
much larger than the typical Overhauser field A/vN 
(with A being the total hyperfine interaction and N be- 
ing the number of nuclei). The timescale tr on which 
RDT describes quantitatively the SE decay well has been 
shown to visibly exceed N/ A. A qualitative difference be- 
tween the SE decay due to interaction with the homonu- 
clear and the heteronuclear bath (containing nuclei with 
distinct Zeeman splittings) is still clearly visible in the 
numerical calculation for 57 w A/vN (corresponding to 
~ 10 mT in large GaAs dots) . Generally, we have found 
that qualitative and even semi-quantitative predictions 
of the RDT are robust down to these fields. The satu- 
ration of the spin echo decay time at low magnetic fields 
had been predicted before using the RDT ) 16 i 17 and our 
new results lead us to expect that this behavior will be 
robust down to B fields corresponding to typical Over- 
hauser field, i.e. we predict the SE decay to be practically 
S-independent between 10 and 100 mT in GaAs quan- 
tum dot with ~ 10 5 nuclei. 

Also, we have found that the SE signal in a homonu- 
clear bath oscillates with the Zeeman frequency of the 
single present nuclear species. This effect is related to 
the well-known ESEEM phenomenon, but to our knowl- 
edge its presence has not been discussed for the spin 
echo signal in case of the isotropic hf coupling. This fea- 
ture might be observed in spin echo experiments on spin 
qubits in quantum dots based on materials having a sin- 
gle nuclear species, e.g. silicon ; 45 ! 46 carbon nanotubes 4 ^ 
or graphene4£ In a heteronuclear bath (e.g. in GaAs) we 
have found that at very low magnetic fields (smaller than 
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the typical Overhauser field) the spin echo signal exhibits 
strong oscillations in which the Zeeman frequencies of all 
the nuclear species are present. 

Finally, we have shown that using a simplified model of 
uniform hf couplings (corresponding to using a box-like 
electron wavefunction) we can recover certain qualitative 
features of the SE signal at timescale exceeding both T 2 * 
and N/A. While the nearly perfect disentanglement of 
electron spin and nuclear bath by the SE sequence was 
discussed previously at high magnetic fields^>21 this re- 
sult of our work suggests that the dynamics of coherence 
recovery in the SE experiment is quite closely related to 
classical dynamics of coupled electron and nuclear spins 
also at much lower magnetic fields. The implications of 
this observation, and also the analysis of accuracy of an- 
alytical theories using the effective hf-mediated interac- 
tion for modeling of the narrowed free induction decay, 
are left for future research. 
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Appendix A: Dependence of the spin echo signal on 
the choice of hyperfine couplings 

In this work we choose the hf couplings Ai randomly 
from a uniform distribution. The realistic distribution 
p(A) of the Ai couplings is related to the shape of the 
envelope wavefunction ^(r) of the electron: 

p(A) = — [ 5(A-A\y(r)\ 2 )d 3 r , (Al) 

where V is the total volume, Vq is the volume of 
the unit cell, and the wavefunction normalization is 
f v | ^ (r*) 1 2 rfr = Vq. For the wavefunction being a two- 
dimensional Gaussian we have p(A) ~ A~ 1 Q(A max — A), 
where O(x) is the Heaviside step function and A max is the 
largest hf coupling (for the nucleus located at the center 
of the wavefunction). p(A) for a more realistic electron 
envelope is given in Ref. [l7|, and it behaves in qualita- 
tively similar way. The key point is that if we are only 
concerned about the most strongly coupled nuclei, with 
Ai not much smaller than A max (which naturally dom- 
inate the decay at short times, and possibly determine 
most of the decay at low B fields), the approximation of 
p{A) by a constant is reasonable. 

The long-time dynamics of the electron spin was pre- 
dicted to be influenced by the shape, specifically the tails, 
of the electron wavefunction ) 15 i 17 i 18 ' 49 ~ — and thus the de- 
tails of the distribution p(A) for small A. These features 



of p{A) are impossible to capture with only N = 20 spins. 
We are, however, not currently concerned about this. 
The recent experiments^, are showing that the low B 
field SE decay occurs on timescale <N/A, on which ac- 
cording to the RDT the wavefunction shape is relatively 
unimportant. Thus the question of low-field accuracy of 
the RDT for an uniform distribution of Ai is well moti- 
vated. 

It is nevertheless prudent to check how our results de- 
pend on the randomly chosen set of Ai. The reason- 
able expectation is that for large N the choice should 
not matter. This is not the case with N < 10, with 
which both the exact and RDT results are showing very 
diverse behavior beyond the short time limit (i.e. i<100 
for N = 10). Unsurprisingly, for such small N the agree- 
ment between the two calculations is also present only 
for these short times. For N — 20 the situation already 
becomes much closer to our expectations. In Fig. [7J we 
present the RDT results for 20 random choices of the set 
of Ai couplings at O = 5.5. One can see that the typical 
decay signal is clearly visible. However, there are still 
sets of hf couplings which give SE signals visibly differ- 
ing from the typical behavior. The two most extreme 
cases are drawn with dashed lines in Fig. [7] The set of 
Ai used in previous calculations corresponds to the dot- 
dashed line. This signal is quite close to the typical one, 
although it does possess characteristic features. These 
features closely correspond to the ones of the exact sig- 
nal shown before in Fig. [2b, showing that at ft = 5.5 
the RDT calculation is very close to the exact one for 
"typical" choice of Ai. 

In Fig. [5] we present the comparison between the RDT 
and exact results for two sets of Ai corresponding to most 
atypical results from Fig. [Jj The values of hf couplings 
used in these calculations are given in Table H together 
with the set of Ai used in most of the other figures in the 
paper. For the signal showing very weak decay (set 1) we 
find quite a good agreement. For the signal exhibiting a 
sharp peak (set 2) the agreement could be considered 
quantitative. 

We have checked that the special behavior in case of set 
2 is due to presence of a large number of very small cou- 
plings in this set. In fact, one can remove more than 10 
smallest Ai from this set without visibly affecting the re- 
sults, while the decay becomes completely different (and 
much weaker) when one of the large couplings, Aig, is 
removed. 

In the case of set 1, the agreement is almost quantita- 
tive (only the slow oscillation seems to be the artifact of 
RDT) As for the reason for the exceptionally weak decay, 
it is possible that after removing the Aj couplings which 
are too small to visibly affect the signal, the remaining 
couplings, especially the ones for different species of nu- 
clei (since hetero-nuclear interactions are crucial for SE 
decay) are too uniform. We note that as shown in Fig. 6, 
the decay for N = 20 with uniform couplings saturates at 
about a half of initial value of the signal. This saturation 
effect becomes weaker when N is increased. 



10 




100 200 300 400 500 600 

t 



FIG. 7: (Color online) RDT calculation of the SE signal 
at fl = 5.5 performed for 20 different sets of hf couplings 
Ai (drawn at random from a uniform distribution, with 
JL = 1 normalization). The signal for the set used in 
previous figures is given by the red dot-dashed line. Two of 
the most "atypical" signals are drawn with dashed lines. 
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FIG. 8: (Color online) The comparison the the RDT calcula- 
tion (dashed lines) and the exact calculation (solid lines) for 
the two sets of hf couplings which give atypical SE signals 
(shown in Fig. [7] by dashed lines) 



It is also interesting to consider how the choice of the 
hf coupling parameters affects the SE signal in the ex- 
perimentally relevant limit of very large N. We consider 
the SE signal S(t), given by 

5(2r) = Tr^e^^e-^^e^^e^] , (A2) 

which is just another way to write Eq. To access 

the limit of large N, we can use the Levy's lemma i 52 ' 53 
It states that for a function f(x) defined on a (N — 1)- 
dimensional hypersphere x £ S l , and satisfying the 
1-Lipshitz condition, the value of the function at a ran- 
domly chosen point is close to the average value of this 
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TABLE I: Three sets of hf couplings Ai: sets 1 and 2 corre- 
spond, respectively, to the most weakly decaying signal and 
to the signal with prominent peak at t ~ 350 shown in Fig. [7] 
while set 3 has been used in all the previous figures. The val- 
ues of uii nuclear Zeeman energies are: Wj = 0.02526 for i < 10, 
cji = 0.0354 for i 6 [11, 16] and ^ = 0.045 for i 6 [17, 20]. 

function (/) with very high probability: 

Prob[|/(f ) - (/) | > e] < cxp (-CNe 2 /L 2 ) (A3) 

where L — sup|V/| is the Lipshitz constant. We use 
the vector (A\,A2, ■ ■ ■ Ak ... , An) of the coupling con- 
stants, with normalization A\ = 1, as a point x on a 
A^-dimensional hypersphere, and the SE signal S(t) as a 
function f(x). To evaluate L, we have to find the deriva- 
tives of 8S(t)/dAk, taking into account that the only 
quantity dependent on Ak in Eq. (|A2|) is the Hamilto- 
nian H . Thus, we need to substitute the equality 

dA k e J Q dA k S tA4j 

= c-^ T f ds(S3 k )(s) 
Jo 

(and its Hermitian-conjugated version) into four places 
in Eq. (|A2[) . In the Equation above (SJfc)(s) = 

e tHs (SJ k )e~ lHs ■ Furthermore, we should take into ac- 
count that the trace of a matrix product Tr[A^B] has all 
properties of the scalar product, so that 

|TrL4tB]| < ^Tt[A^A}Tt[BW] = ||i||||B||, (A5) 
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where 
dS(t) 



dA k 
Therefore 



^/Tr[^4t^4], As a result, we obtain 



•14||.S':||/ dfl||(SJ fc )(*)|| = V2 / ds||(SJ*)(«) 
o Jo 

(A6) 



dS{t) 



OA, 



< C\t 



(A7) 



with some constant C\ independent of iV, and for the 
Lipshitz constant we obtain 



L < CitVN. 



(A8) 



Thus, for t ~ 1 (in our dimensionless units, where 
Sfe^fc = 1)> the specific choice of the hf coupling pa- 
rameters does not matter. 

However, it is obvious that our rigorous estimate 
Eq. (|A6|) . which is based on straightforward use of the 
Cauchy inequality, is too crude. Heuristically, we expect 
that the correlation between the central spin and a sin- 
gle bath spin, which we crudely estimated from above as 
0(1), is actually of order of 0(1/N). Then, the range 
of times where the choice of the hf parameters does not 
matter, is extended to t ~ N. 



As discussed in Sec. IIV1 the Hamiltonian couples only 
pairs of |±,j, m) states: 



-iHr 



T \+J,m) = a jm \+,j,m) + b jm \-,j,m + l) ,(B3) 
\-,j,m) = c jm \-,j,m) +d jm \+,j,m- 1) (B4) 



and the coefficients aj m and bj m are given by 



m 



N+ 



srn^r (B5) 



V+ AJ + 

■ -iE + t jm ■ jm 

-te m — ^r— sin — - — r 
N + 2 



(B6) 



and Cjm and dj m are given by analogous expressions with 
superscripts + replaced by — . E m , X ± ,Z ± and N ± are 
given by 



E 
X 
Z 



± 

m 

± 
jm 

± 



[(2m±l)u]-A/2}/2 , (B7) 

Ay/j(J + l)-m(m±l) , (B8) 

±[Q-u) + A(m±l/2)} , (B9) 

(Xf m f + {Zk? . (BIO) 



Appendix B: Spin echo in a model with uniform hf 
couplings 

We rewrite Eq. (01 using the basis of |±,j, m) states 
(with electron states \±x) — (|+) ± |— ))/\/2) and with 

T = t/2 



x m 4eedj 



J rn=-j 



(+x,j,m\e intT a x e inT ^-e~* nT a x e~ inT \+x,j,m)(l}l) 

where the sum over j is from to N/2 (1/2 to N/2) for 
even (odd) number of nuclei N, and the degeneracies of 
states with given j are given by^ 



D, 



N\ 



2j + 1 



(JV/2-j)!(JV/2+j)!iV/2 + i + l 



(B2) 



Plugging these into Eq. (|B1[) and into an analogous 
expression for (S y (t)) (in which the middle a x operator 
is replaced by a v ) we arrive at the formula for the deco- 
herence function 



N/2 3 



W SE (t) = J2 J2 ^{\ a J^c ]m \ 2 + a ]m c* m _ 1 \d jm \ 2 + 

j=0 m=-j 



\bj m ?) 



(Bll) 



In the case of heteronuclear bath with Nj nuclear 
species one has to introduce Nj sets of basis states 
\j a ,m a ) in Eq. (|B1[) . The Hamiltonian is then coupling 
the whole subspaces of fixed j a , and one has to consider 
Hamiltonian matrices of dimension 2rj a (2j Q + 1) and 
evaluate the evolution operators numerically. 
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